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Recently, interesting nonmonotonic time evolution has been pointed out in the experiments 
by Jonsson et al. and Jonsson et al. and also in the numerical simulation by Takayama and 
Hukushima, where the magnetic susceptibility does not monotonically relax to the equi- 
librium value, but moves to the opposite side. We study the mechanism of this puzzling 
nonmonotonic dynamical property in a frustrated Ising model in which the equilibrium cor- 
relation exhibits nonmonotonic temperature dependence (reentrant type). We study the time 
evolution of the spin correlation function after a sudden change of temperature. We find that 
the value of the correlation function shows nonmonotonic relaxation, and we analyze the 
mechanisms of the nonmonotonicity. We also point out that competition between different 
configurations generally causes nonmonotonic relaxation. 

KEYWORDS: frustration, decoration bond, reentrant phase transition, relaxation, Ising spin 
system 

Recently, Takayama and Hukushima pointed out an interesting dynamical property of the 
magnetization of Ising spin glasses in the magnetic field after the halt of the field cooling. 1 
They studied the dynamics of magnetization at various temperatures and in various field pro- 
tocols. That is, they decreased the temperature at several speeds in a finite magnetic field, and 
observed the change in magnetization over time. In particular, they found that the magneti- 
zation after the halt of field cooling beyond the spin glass transition temperature sometimes 
shows nonmonotonic relaxation. In some cases, the magnetization at the end of the field cool- 
ing process was smaller than the equilibrium value. However, the magnetization still decreased 
which was the opposite tendency to the equilibrium value. The observation is consistent with 
the experiments of the field-cooled protocol in Fen.sMno.sTiOs 2 and Feo.55Mno.4sTi03, 3 which 
are good model systems for a short-range interacting Ising spin glass. This type of nonmono- 
tonic dynamics in the relaxation toward the equilibrium state is interesting, and in this Letter 
we study possible mechanisms using simplified models. 

In frustrated spin systems, the correlation often shows nonmonotonic dependence on the 
temperature. 4,5 This nonmonotonic dependence causes reentrant phase transitions, 5-7 and also 
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can be an origin of temperature chaos in spin glasses. ' In a simple ferromagnetic model, a 
phase transition occurs between the high-temperature paramagnetic phase and the ferromag- 
netic phase. The ferromagnetic state is energetically favorable and realized at low temperature. 
This situation gives a standard order-disorder phase transition. In contrast, in frustrated sys- 
tems, different types of ordered states are often nearly degenerate. In such cases, different 
order phases are realized at different temperatures. For example, when the temperature de- 
creases from a high temperature, the system may show successive phase transitions, i.e., from 
the paramagnetic to an ordered phase (e.g., with an antiferromagnetic order) and then to 
another ordered phase (e.g., with a ferromagnetic order). Generally, this type of successive 
phase transition is called a reentrant phase transition, although originally "reentrant" means 
the phase transitions from paramagnetic phase to an ordered phase and then to the param- 
agnetic phase again. The nature of reentrant phase transitions has been studied exactly using 
the two-dimensional Ising model. 6 This nonmonotonic dependence is not only important in 
the phase transition but also plays an important role in the temperature dependence of the 
short-range correlations. In frustrated systems such as spin glasses, various frustrated local 
interaction configurations are realized. In such systems, we expect the local correlations to 
show nonmonotonic temperature dependence, and we expect the ordered structure to change 
globally with temperature, which can be a possible mechanism of the temperature chaos and 
the rejuvenation phenomena. 8 



Fig. 1. (a) A frustrated lattice. The dotted line denotes Jaf and the solid lines denote Jp. The 
effective coupling between o\ and a 2 is given by K^(T). (b) A frustrated lattice in which the 
structure of Fig. 1(a) is used n times. The effective coupling in this structure is K^(T). 

In this Letter, we study the dynamics of the correlation function in a frustrated Ising 
model 



which shows reentrant-type temperature dependence of the correlation function. Because of 
the freedom of the local gauge transformation {&i — > —ai,Jij — > —Jij}, we can have many 
configurations of interactions with the same property of frustration. We call this degeneracy 
"Mattis degeneracy". In the case of H = 0, all the models in the Mattis degeneracy have 




(a) 
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the same thermodynamic properties. On the other hand, the magnetic response, of course, 
depends on the configuration {Jij}. Thus, we mainly study the correlation function of spins, 
which is essentially the same (apart from the sign) in the Mattis degenerate systems. 

One of the simplest models of this type is depicted in Fig. 1(a). The Hamiltonian is given 

by 

H = Jaf<7i02 - Jf0"3(oi + 0"2), (2) 

where 2Jp > Jaf > 0. The effective coupling K e s(T) between the spins o\ and 02 at temper- 
ature T is defined by 



<7 3 =±1 

and effective coupling between o\ and 02 is 



A(T)e K ^ {T)(Tl(T2 



■ , H — log ( cosh f - — ^ ) , 
k B T 2 6 V \^bT))' 



(3) 



(4) 



where A(T) is an analytic function of T, and the equilibrium correlation (0102) is given by 
tanh-Keff . In Fig. 2, we plot the equilibrium correlation (0102) as a function of the temperature 
for the parameters Jf = 1 and Jaf = 0.5. The bold solid curve denotes the case of the model 
given by eq. (2). Hereafter, we take Jf as the unit of energy. The temperature is also scaled by 
Jf. In Fig. 2, we can observe nonmonotonic temperature dependence. The correlation (0102) 




Fig. 2. (Color online) Temperature dependence of the correlation function (0102) for n =1, 2, 4, 6, 
8, and 10. The bold solid curve denotes the data for n = 1. and the bold dotted line denotes the 
data for n = 10. The intermediate lines show n=2, 4, 6, and 8. 



at the higher-temperature side is small. To increase this amplitude, we provide a multiplied 
decoration bond with n intermediate spins (the open circles depicted in Fig. 1(b)). Therefore, 
the total number of spins is N = n + 2. The Hamiltonian in this case is given by 



H {n) = nJ AF 0i0 2 



n+2 
fc=3 



(5) 
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and the effective coupling between o\ and 02 is given by 

^m-^cn— ^ + =^(«h(g)). (6) 

The correlation functions for various values of n are also plotted in Fig. 2. 
Now we study the dynamics of the correlation function 

C ( t ) = E P ({^},t)cna 2 , (7) 

{<Ti=±l} 

where P({ai},t) is the distribution function at time t. We adopt the Glauber-type kinetic 

Ising model 10 for the time evolution 

dP(ai,--- ,Q-j,-- - ,<7N,t) 
dt 

= ~y^P(o-i,--- ,CTj,--- ,PN,t)Wai^-<Ti + E- P ^ CJl ''" > _(T *>'" J^iV,*)^-^^, (8) 

with the transition probability per unit time 

Peq(oi) • • • j — Cij • • • , CTat) 



lb 



• • • 5 f at) ' 



where 



-Peq(0"l, • • • , <Ti 



p -/9W({o-i}) 



(9) 



(10) 



It is convenient to use the vector P(t) consisting of the probabilities of the states 

/ P(++,-- - , + + +,t) \ 
P(++,--- , + + - t) 



P(t) 



The dynamics is expressed by 



V n- 



P(t + At) = CP(t), 



(11) 



where £ is a 2^ x 2^ matrix with matrix elements at the small limit of At. 



Cij = -j^wj^iAt for i ^ j 



(12) 



(13) 



and 



Ca 1 ^ Cp. (14) 

Hereafter, we adopt eq. (12) as an approximate discrete form of eq. (8). Here, we suddenly 
change the temperature from T\ = 10 to T 2 = 1 and study the subsequent dynamics by 
iterating the time evolution operator C with At = 1. We call an update by this procedure a 
"Monte Carlo step (MCS)". The initial probability distribution is set to be the equilibrium 
one at T = T± . 
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Fig. 3. (Color online) Time (Monte Carlo step) dependence of the correlation function for n = 2, 4, 
6, 8, and 10. The types of line mean the same as those in Fig. 2. 

In Fig. 3, we depict the time evolution of the correlation C(t). We find that C(t) first 
decreases, the amplitude increases on the negative side, it then becomes positive, and fi- 
nally it reaches the equilibrium value at T = Ti- This observation indicates that even if the 
temperature is changed suddenly to T2, the correlation function does not necessarily relax 
directly toward its new equilibrium value, but it can show nonmonotonic relaxation. We may 
understand that the "effective temperature" of the system decreases gradually even when 
the temperature of the thermal bath is changed suddenly. The correlation function, as well 
as other quantities, shows a similar dependence to its temperature dependence. That is, if 
a quantity shows a nonmonotonic temperature dependence in equilibrium, it tends to show 
nonmonotonic relaxation. 

To understand this nonmonotonic behavior, we analyze the time evolution from the view- 
point of the eigenvalue problem of the time-evolution operator C. 11 Let u k and X k be the /cth 
right eigenvector of C and the eigenvalue, respectively (k = 1, • • ■ , 2^) 



Here, we assume that 1 = Ai > A2 > A3, • • • , > A 2 jv. The state u\ gives the equilibrium state 
at T = T%. The initial state (the equilibrium state at T = Ti) is expanded using {u k } as 



£u k = X k u k . 



(15) 




(16) 



k=2 



After t time evolutions by jC, the state evolves to 




(17) 



k=2 



The contribution of the kth mode to the correlation function is 




(18) 



{<Ti=±l} 
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The sign of Ck can be either positive or negative, and the contribution of each mode Ck{t) 
causes a decrease or increase in C(t) corresponding to the sign of Ck S{ CTi =±i} u k({ cr i}) cr l cr 2- 
The relaxation time of each mode is different and C(t) can be nonmonotonic. Here we demon- 
strate this situation taking the simplest case of n = 1. 




-0.4 - 

o 

I , , , , , , , , , 
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Fig. 4. (Color online) (a) Eigenvalues Xk for the case of n = 1 (solid circles) of the time evolution 
operator at T = T2(= 1), and the coefficients Ck (open circles) for the equilibrium state at T — 
Ti(= 10). Only the modes of k = 4 and 7 give contributions, (b) Time evolution of C(t) (bold 
curve) with the contributions from the 4th mode (solid diamonds) and 7th mode (open diamonds). 
Open circles correspond to the sum of c^. 



In the case of n = 1, there are 8 modes. The eigenvalues and the coefficients c& are 
plotted in Fig. 4(a). We take the initial state to be the equilibrium state at T±, and the 
contribution from some modes is zero because of the symmetry. For example, modes that are 
antisymmetric in exchanging sites 1 and 2 do not contribute. In the present case, only two 
modes k = 4 and 7 contribute. The changes in C±{t) and CV(t) are depicted in Fig. 4(b). The 
contribution from the 7th mode with a short relaxation time (open diamonds) is positive and 
it relaxes rapidly. This causes the decrease in C(t) in the early stages. On the other hand, the 
contribution of the slow-relaxing mode (A; = 4) is negative and it causes an increase in C(t) 
to the equilibrium value. For larger n, we found similar behavior (not shown). 

It is interesting to note the characteristics of the relaxation modes. In Fig. 5, we plot the 
probability distribution in the phase space of the equilibrium states P eq (T) at T = 1, 2, and 10. 

The configurations 1, 2,3,4, • • • ,8 denote (<t 3 , a 2 , oi) = (+ + +), (+ + — ), (H h), (H ), 

• • • ( ), respectively. In configurations 1, 4, 5, and 8, the spins a\ and 02 have the same 

sign and thus there is a ferromagnetic correlation between a\ and and the configurations 
2, 3, 6, and 7 have an antiferromagnetic correlation. We also plot the eigenvectors and 
U7 at T = T 2 {= 1). There, we find that the fast relaxation mode 117 has large amplitudes 
at configurations 4 and 5, and the reduction of this mode causes the change from P eq (10) 
to P cq (2), where the probabilities of the energetically unfavorable configurations 4 and 5 
are reduced, and the antiferromagnetic correlation between o\ and a 2 slightly increases. The 
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Fig. 5. (Color online) Probability distributions of the equilibrium states at T = 1 (open circles), 2 
(squares), and 10 (closed circles), and eigenrelaxation modes U4 (crosses) and uj (pluses). Lines 
between points are drawn to highlight the trends. The configurations 1,2,3,4, •• • ,8 denote the 
states (0-3, (72, 01) = (+ + +), (+ + -), (H h), (H ) • • • , ( ), respectively. 

slow relaxation mode 1x4 corresponds to the difference between P eq (2) and P eq (l) where 
the ferromagnetic correlation between o\ and 02 increases. Here we find that the relaxation 
modes reflect the temperature dependence of the equilibrium state. This fact corresponds to 
the gradual cooling of the system. 

Next, we study the antiferromagnetically coupled spin cluster in a uniform magnetic field 

Ti = J(aia 2 + 02O3 + 0-30-4 + 0-40-1) - H{a\ + a 2 + 03 + 04), (19) 

where J and H are competing. In other words, this system has "frustration" between the 
magnetic coupling and the magnetic field. In Fig. 6(a), we depict the temperature dependence 
of the magnetization M = o\ + 02 + 0-3 + 0-4 in the magnetic field H = 0.1J. The magnetization 
shows nonmonotonic behavior as a function of temperature. The temperature is scaled by J. 
The time evolution of M after sudden quenching of the temperature from T = 10 to 0.1 is 
plotted in Fig. 6(b). Here we again find the same nonmonotonic behavior. 




(a) 




(b) 



Fig. 6. (a) Temperature dependence of magnetization M of the spin antiferromagnetically coupled 
spin cluster (eq. 19). (b) Time-evolution of M after sudden quenching of the temperature from 
T = 10 to 0.1. 
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We also refer to nonmonotonic behavior in macroscopic models. If the system has a 
metastable state that is separated from the equilibrium state by a free energy barrier, as 
depicted in Fig. 7, the system can show nonmonotonic behavior. For example, if the initial 
state is point B, it relaxes to the metastable point in a short time, and later it relaxes to the 
equilibrium state via a nucleation process. This is characteristically different from the simple 
relaxation starting from point A. Thus, this type of difference in the relaxation has been used 
to detect a metastable state. 12 In the present macroscopic case, there are also two competing 
structures representing the metastable state and the equilibrium state. 




stable 



Fig. 7. Schematic free energy structure with a metastable state. 



In the present study, we only investigated the cases with sudden changes in the tempera- 
ture. It is also interesting to change the temperature with a finite sweep velocity. If we sweep 
the temperature sufficiently slowly, a quasi-static state should be realized. The sweep velocity 
dependence of the dynamics will be studied in the near future. 

Using a lattice constructed by the decoration bond, we can study the reentrant phase 
transition. 8,9 If we change the temperature suddenly from the paramagnetic region to the 
ferromagnetic region, we obtain very slow relaxation due to an entropy-induced screening 
effect. 9 Here, we find that the antiferromagnetic correlation appears for only a very short 
time and the nonmonotonicity does not develop. To achieve nonmonotonic behavior on a 
macroscopic scale, we have to construct a system where the present mechanics takes place 
on a coarse-grained scale. Here, the present spin Oi should represent a coarse-grained local 
magnetization. In spin glasses, the frustrated configuration remains in each step of the coarse- 
grain magnetization process, which may be called a "hierarchical structure". Furthermore, 
in real spin glasses, the spin directions in the ordered configuration are spatially random, 
and we need to average these over all the possible configurations. Thus, the temperature 
dependence of the uniform magnetization is averaged out, although the local spin correlations 
may show various peculiar temperature dependences. However, the hierarchical frustrated 
structure causes peculiar phenomena in the time dependence. As was proposed by Takayama 
and Hukushima, the magnetization first decreases while the short-range spin glass correlation 
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still increases, whereas it finally increases after the correlation reaches the saturated value 
where the local magnetization behaves independently. This complexity causes the peculiar 
spin glass properties, but the present mechanism may give the basic local mechanism of these 
phenomena. 
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